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ABSTRACT 

We investigate the effects of intermittent strong jets from an Active Galactic Nuclei 
(AGN) of a massive galaxy in the core of a cool core galaxy cluster, on the dynamics 
and metal distribution of the intra-cluster medium (ICM). We use a simple model 
for the metal distribution within the host galaxy which includes metal injection via 
star formation. We carry out 2D axisymmetric hydrodynamic simulations of these 
systems. After having established a cooling flow, two light jets are injected in opposite 
directions with a range of (intermittent) active phases. We follow the time evolution of 
the system from the jets' active phases up to 3 Gyr. The general metallicity evolution 
for all our simulations is very similar on large-scales. The convective flows driven by the 
jets advect gas and metals from the central galaxy to distances beyond 1.5 Mpc within 
the cluster. Intermittent jets are able to distribute metals to grater radii. The metal 
injection has effects on the general metal abundances present in the ICM, the faster 
the metal replenishment time the higher the metal abundance. Although metallicity 
gradients in the very central regions of cool core clusters are likely to be shaped by 
less-energetic phenomena, we find evidence in our simulations for metallicity gradients 
similar to those observed out to ~400kpc radius. The metal distribution details of 
the central galaxy have no effects on the dynamical evolution of the ICM metals. 

Key words: hydrodynamics - galaxies: jets - methods: numerical - X-rays: galaxies: 
clusters. 



1 INTRODUCTION 

The abundance and distribution of heavy elements (met- 
als) in the intra-cluster medium (ICM) holds information 
about the history of galaxy clusters (clusters from now on). 
Metals formed as a result of star formation in the cluster 
galaxies must be transported into the ICM and a variety of 
mechanisms are possible, for example galactic-scale winds, 
or outflows driven by Active Galactic Nuclei (AGN). X-ray 
spectroscopic observations provide a direct prove of metal- 
licity gradients in the ICM. Many clusters show an approx- 
imately flat metallicity distribution, of about a third of the 
solar (.Z(7)), out to megaparsec d istances from the cluster 
core jPe Grandi fc MolendillioOll ). Moreover, the metallic- 
ity profiles of cool core intermediate redshift (~ 0.1 -0.35) 
clusters, particularly those that show no evidence of shocks 
associated with mergers, have a n abundance excess of the or- 
der 0.45-0.55 Z© at their cores (|De Grandi fc Molendl200ll : 



* E-mail: mh475@mrao.cam.ac.uk 



iBaldi et all l2007l : iLeccardi fc Molendil [20oi ). High resolu- 
tion Chandra X-ray observations of the core of Hydra A, 
show that metals in the ce ntral re gions were formed by dif- 
ferent stellar populations (jPavid e t al. 2001), with a metal 
replenishment time of the order of 60 M© yr~^ over a period 
of 1 Gyr. Simil ar results are found for other clusters as well 
(|Molendill20"0i l. 

Over recent years there has been a growing consensus 
that radio galaxies h ave a significant impact on the clus- 
ter environment (see iMcNamara fc NulsenI |2007| for a re- 
cent review). Jets from radio- loud AGN provide an ideal 
mechanism for coupling the energy output from the AGN 
to the surroundirig gas, both via he ating by the bow shock 
iBinnev fc Taboj|l995l : lKaiser fc A lexander IQQ j Reynold s, 
Heinz fc Begelman 2001: lAlexmden ,2002 : ,Kraus3 I2003I '). 
but also by driving gas from the core during both the ac- 
tive phase of the jet (e.g., Alexander 2002), but also as 
the buoyant remnant - bubbles - o f the radio-source co- 
coons rise though the cluster gas (e.g.JChurazov et "al]l200ll : 
iBlanton et allbOOll : [Reynolds etaL^l2001^ . The effects can 
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be very long lasting as a large-sc ale convective flow can be es - 
tablished by the rising bubbles l|Basson fc Alexandeij[20o3 ) . 
These flows have important effe cts on the entropy and en- 
ergy of the central ICM (see e.g. jKaiser et al1l2005l and ref- 
erences there in), and can transport metals from the core 
to ra dii of order a megaparsec within the cluster ijRenzinil 
I2OO4I : Heath, Krause & Alexander 2007). Moreover, the bub- 
bles are i nferred to produce h ydrodynamic turbulence in 
the ICM (iRebusco et al.ll2006l ). These processes may well 
be linked to the cluster-wide feedback with infall of coolini 



gas triggering the AGN activity (e.g., Binnev fc Taboi]ll99_. 
iKaiser fc Binnevl[200^ : IShabala et al.ll2008l ). The observed 
lack of evolution of the ICM metallicity out to redshift of 
order 1 JCariUi et al l 200 ij : iMushotzkv fc Loewensteinll 19971 : 
iLeccardi fc Molendil 20081 ). suggests that metal advection 



took place at or before a redshift of unity. Furthermore, the 
metallicity gradients in the ICM also suggest that the met- 
als were expelled from the dominant central cluster galaxy 
rather tha n from ram-pr essure stripping of cluster galaxies 
m general (iRenzinilbOOi ). Supernovae do not provide suffi- 
cient energy for that proce ss and an attractiv e alternative 
are powerful radio galaxies (|Heath et al.ll2007l ). 

In this paper we examine the advection of metals driven 
by radio- loud AGN. We investigate numerically the effects 
that the powerful jets and continuous metal formation, of a 
massive galaxy in the core of a cool core cluster, have on the 
hydrodynamics and chemistry of the ICM. We explore the 
consequences of varying the jet power, the galaxy metal dis- 
tribution profile and the metal formation rates. Our imple- 
mented cluster and galaxy models are described in Section 2 
along with our numerical setup. We describe and discuss our 
results in Section 3 and give our co nclusions in Section 4. 
In this paper we extend the work of iHeath et al.l (j2007l ). In 
particular, we use an improved numerical scheme (see be- 
low) which enables us to model intermittent jets as well as 
injection of metals from star formation. 



2 THE SIMULATIONS 

2.1 Numerical setup and initial conditions 

We so lve the equations o f hydrodynamics using the Flash 
code ijFrvxell et all l2000l l assuming axisymmetry in a 2D 
sphere. Given this geometry we us e spherical polar co ordi- 
nates (r,6'). Flash is a well tested ()C alder et al.|[2003 ) par- 
allel modular hydrodynamics code with an adaptive mesh. 
It solves the Riemann problem using a Piecewise-Parabolic 
Method able to work well with contact discontinuities, mak- 
ing it well-suited to our current problem. 

To implement the galaxy cluster we use a computational 
domain with ranges r £ [10,3000] kpc, with outflow and hy- 
drostatic boundary conditions at the inner and outer radii 
respectively, and 6 G [0,7r] rad both with reflective boundary 
conditions to establish the rotational symmetry of the sys- 
tem. Our initial grid setup consisted of 512x16 zones which 
were adaptively refined by a half, attaining a maximum res- 
olution of 2.92 kpc in the radial direction and 1.27 kpc in the 
angular one. Each simulation took 5 hours (on average) to 
run on 20 processors of the CamGritfl cluster of the Univer- 



sity of Cambridge, and the production runs were made at 
the Darwir0 cluster of the University. 

The initial conditions of our ICM gas were simulated as 
a monoatomic hydrogen plasma with an ideal gas equation 
of state (7 — 5/3), an initial constant temperature of 10^ K 
and a King density profile ((Kine 1972. '): 



Pair) 



Po 



(1 -f (r/ao)2)3/3/2 



(1) 



We set the cluster's core density, po, and radius, ao, to be 
10~^^ kgm~^ and 100 kpc respectively, and take (3 = 2/3. For 
this initial gas distribution to be in hydrostatic equilibrium, 
the gravitational acceleration of the dark matter potential 
well must satisfy: 



9r 



ZPkBT 



p,mHal (1 + {r/aoY 



(2) 



(3) 



where fj,= 1/2 and ttih is the proton mass. We neglect the 
self-gravity of the cluster gas and take this as a fixed back- 
ground gravitational field. 

We include an energy loss term, dE/dt, to account for 
the ICM radiative cooling due to Bremsstrahlung: 

dE/dt = -1.722 X 10-« (-^)' ^^W/m-. 

Our simulations consist of three phases: (1) The cool- 
ing fiow phase, in which we let the initial ICM to evolve 
subject to Bremsstrahlung radiation and gravity only, until 
it reaches a quasi-steady state. The system then develops 
a nearly spherically symmetric inflow. (2) The radio galaxy 
active phase(s), when we use the 9^ [0.5,7r— 0.5] rad region 
of the inner radial boundary (with r= 10 kpc) to inject two 
light jets into the cluster with a velocity of a third of the 
speed of light (Lorentz factor = 1.048). We set the density of 
the jets to be 10~^ times smaller than the initial central ICM 
density, and the pressure of the jets was adjusted to that of 
the cooled central ICM of 2.2 x 10~^'^Nn i~^ This simulates 
a powerful FR II (jPanaroff fc Rilevlll974h radio source with 
a jet power of 8.8 x 10'^^ Watts. (3) The remnant phase, dur- 
ing which the ICM, the metals (see the following Section) 
and the convective flows left by the jets evolve passively, un- 
der the influence of gravity and Bremsstrahlung radiation, 
for 3 Gyr. Our different simulations and parameters are sum- 
marised in Table 1. 



2.2 The tracer fields, conservation and injection 

We implement the distribution and evolution of metals using 
Flash's tracer flelds (fields from now on). The metallicity, 
Z, is the ratio of the metal density to the ICM gas density, 
Z{r) = pm.{r)/pg{r), and the field is defined so that: 



f{r) = NZ{r), 



(4) 



where A*' is a constant factor required by Flash to keep 
/ £ [0,1]. Flash uses two flelds: /, representing the metal 
distribution, and /, the non-metal gas fraction. Then, flelds 
are conserved in the sense that / 4- / = 1 in every zone 
of the simulation grid. Thus, normalisation parameters are 
required to fulfll this condition. 



. escience. cam. ac. uk / proj ects / camgrid / 
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Table 1. Simulations and parameters. The names of the simulations contain the information 
about the time between the jet outburst (interludes, I), the metal injection replacement time 
(R) and a letter P when the central galaxy metal distribution is implemented with a Plumer 
profile. 



Simulation 


Number 


Duration of 


Time between 


Injection re- 


Metal 


name 


of jet 


each outburst 


outbursts 


placement time 


distribution 




episodes 


[Myr] 


[Myr] 


U [Myr] 


profile 


lO-RInf 


1 


30 





oo 


flat 


lO-RlOOO 


1 


30 





1000 


flat 


I0-R500 


1 


30 





500 


flat 


IlO-RInf 


3 


10 


10 


oo 


flat 


IlOO-RInf 


3 


10 


100 


oo 


flat 


IIOO-RIOOO 


3 


10 


100 


1000 


flat 


I100-R500 


3 


10 


100 


500 


flat 


IlOO-RInf-P 


3 


10 


100 


oo 


Plumer 


IIOO-RIOOO-P 


3 


10 


100 


1000 


Plumer 


I0-R500-P 


1 


30 





500 


Plumer 


I100-R500-P 


3 


10 


100 


500 


Plumer 



We use two galaxy metal density distributions: a flat 
distribution: 



for r £ [10, 50]kpc; 
elsewhere. 



and a Plumer- like distribution (see Table 1): 



Pmp{r) 



(l-H(r/r,)2)5/2' 



(5) 



(6) 



where rc in the denominator of equation (|6]) is the core ra- 
dius of the galaxy metal distribution, which we set to 15 kpc. 
The metal density, p^p , in equation (|6} was adjusted to be 
equal to po (equation [ij making the metal density, Pmp{r), 
a reasonable and arbitrary fraction of the ICM gas. On the 
other hand, pgj,, in equation ([Sjl, was adjusted so that the 
initial metal mass in the cluster is the same for all our sim- 
ulations. 

Metals are injected via star formation which we con- 
sider only in the central dominant cluster galaxy. To simu- 
late this process (from phase 2 of the simulation onwards) we 
add metals to the ICM gas, with a radial dependence follow- 
ing the initial metal distribution -either Po/('") '^^ pmpir) 
(equations [S] or respectively) - with a metal replacement 
time, tc, so that: 



dpm{r,t) _ pm{r) 
dt tc ' 



dpgjr, t) ^ pm(r) 
dt tc 



(7) 



Values for tc (see Table 1), are based on p revious work which 
consider a constant star formation rate I David et al.l l200ll : 
Gaibler, Krause & Camenzind 2005. From equations ([7]) the 
computational Euler step for these densities is given by: 



(r,t)= pUr,t)+prn{r)(dt/tc), 
p''+\r,t) = p"g{r,t) + p^{r){dt/tc), 



(8) 



where n + 1 denotes the densities at the current timestep, dt, 
and n, the densities at the previous timestep. The evolution 
of the metallicity is hence given by equations and 



+ Pm{r) (dt/tc) 



^ ' P",+' " p^+p^{r) (dt/tc)' 
and thus the evolution of the fields follows so that: 



(9) 



Npmjr) (dt/tc) 
Pa + PAr) [dt/tc) J + Pm{r) (dt/tc) 



(10) 



where /"'''^ and /" represent the fields for the current and 
the previous timesteps respectively. 

We tested our implementation and conservation of the 
fields by following the evolution of an annulus sector with a 
constant field distribution, during a cooling fiow phase (see 
Section 2). We found a loss of fields due to advection of 
2 parts in 10^ over a Hubble time. The inner computational 
boundary has outflow boundary conditions (except for the 
jets' nozzles during the active phases) and thus fields (met- 
als) can be lost only if they fiow through this boundary. 
A similar test for our implementation of fields injection re- 
sulted in a metal loss of not more than 1 %, over a Hubble 
time - since our simulations last for 3 Gyr, we confidently 
neglect these numerical effects. 



3 RESULTS AND DISCUSSION 

The basic dynamical evolution observed in the simulations 
is as follows. During phase 1 (see Section 2) radiative cooling 
leads to the formation of a cooling flow. We run the simula- 
tion until the cooling flow reaches a steady state and devel- 
ops a core with cool inflowing gas and a radial gradient in 
the velocity - this gradient is very steep near to the core and 
flat away from it. The initial hydrostatic equilibrium condi- 
tions prevail away from the core in agreement with canonical 
models (e.g., see Fabian 1994). During phase 2, the initial 
hydrodynamical structure of powerful jets develops. The jets 
drive a supersonic bow shock from the centre of the galaxy 
through the ICM, but at this stage the flelds are not removed 
from the core. Behind the shock, the under-dense jets form 
a low density cocoon bounded by a contact surface. The jets 
push aside the central cool gas and metals in their way. The 
cocoon g rows roughly self-sim i larly in ac cordance with th e 
models of lKaiser fc Alexandeil l|l997h and lAlexande^ l|2006l ). 
The left panels in Figure [T] show the distribution of the log- 
arithms of the density (top) and the fields (bottom), at the 
end of the active phase of run lO-RInf (30 Myr). The jets, 
the cocoon and ambient medium are clearly visible. 
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Figure 1. Plots (xy plane) of the logarithms of the ICM gas density (top) and the tracer fields (bottom). Left panel shows run lO-RInf 
(see Table 1) at the end of the jets' active phase. The jets and the cocoon are visible. The middle and right panels show runs lO-RInf 
and I100-R500 respectively, at the end of the simulation (3Gyr). 



Once the cycle of the jets is finished, the system has 
developed several gas layers with different densities and ve- 
locities. The expansion of the cocoon decelerates gradually. 
Hydrodynamic instabilities grow and disrupt the cocoon 
producing turbulent flows in the central ICM. The rem- 
nant cocoon forms two hot buoyant bubbles that propagate 
away from the core, in agreement with previous 3D hy- 
drodynamical simulations (e.g., see IChurazov et al.1 12OOII : 
iBasson fc Alexander! [2003 ). 

The fields mimic the motion of the gas fiows closely 
and mixing of gas with low and high metal concentrations 
occurs, particularly close to the most prominent turbulent 
vortices. The later effects are more important when inter- 
mittent jets are present (see Table 1). The convective flows 
drag the metals very efficiently - ~ 95% of the initial galaxy 
metals in general - to distances up to 2 Mpc (see Figures [T] 
and [2]). The bow shock continues propagating, decelerates 
gradually and reaches the outer computational boundary 



(3 Mpc) at the end of our simulations (3 Gyr) . The middle 
and right panels in Figure [l]show the distribution of the log- 
arithms of the density (top) and the fields (bottom), of these 
convective fiows, at the end of the simulation for Runs 10- 
Rlnf and I100-R500. 



3.1 The time evolution of the metal distribution 

Figure [2] shows the (logarithm) azimuthal average of the 
relative metallicity of the ICM gas density, as functions of 
radial distance from the cluster's core at different epochs. 
The relative metallicity in a given shell of thickness dr is 
calculated as: 

Z(r) _ 1 J pr4r',9)r"'sin{e)dedr' 
Zo ^ Zo f pg{r',ey^sin{9)d9dr' ' 

where Zq is an arbitrary normalisation factor which, in Sec- 
tion [3]2] below, we use to compare our simulated metallicity 
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Figure 2. Time evolution of the mean relative metallicities as functions of the radius from cluster centre. Each point shows the logarithm 
of the azimuthal integration (spherical shells) of the ratio of the metal density over the ICM gas density (see equation [Tl]l . i.e., the mean 
metallicities are normalised in order to be compared with the observations (see Section |3.2| I. Column 1, comparison of simulations of 
galaxies with no metal injection and different jets' duty cycles. Column 2, comparison of Runs with a single 30 Myr outburst and different 
matter injection rates. Column 3, comparison of simulations with three 10 Myr outbursts (and two 100 Myr interludes) and different 
mater injection rates (see Table 1). 
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Figure 3. A zoom into the central 300 kpc of Figure [2] 



gradie nts with the observations from iLeccardi fc Molendil 
(20081 ). The left most column in Figure [2] (hereafter 2L) 
shows Runs without any metal injection and with different 
jet active phases. The middle column (2M) shows Runs with 
a single 30 Myr outburst and different metal distributions 
and injection rates. The right most column (2R) displays 



Runs with the same episodic jets, but different injection 
rates (see Table 1). Furthermore, Figures [3] and [J] present 
plots with a similar format, arranged by columns (and we 
address them with the same short notation for the left mid- 
dle and right columns). Figure [3] shows details of the central 
region of all the plots in Figure [2] while Figure |4] compares 



© 2007 RAS, MNRAS OOO.nifTOl 



Dynamical and chemical effects of FR II radio sources on the intra-cluster medium 7 



No Injection, Episodes 100 iVIyr 

lOO-RInf 
100-Rlnf-P 



500 



2000 



1000 1500 

Radius ri<pcl 
No injection, Episodes 1000 IVIyr 




500 



2000 



1000 1500 
Radius [l<pcl 

No Injection, Episodes 3000 Myr 




1000 1500 
Radius [l<pc] 




Singie Jets, Inject. 100 Myr 

I0-R500 
I0-R500-P 



500 



2000 



1000 1500 
Radius fkpc] 
Single Jets, Inject. 1 000 Myr 




500 



2000 



1000 1500 
Radius [kpc] 
Single Jets, Inject. 3000 Myr 




1000 1500 
Radius [kpc] 



Episodic Jets, Inject. 100 Myr 

noo-Riooo 

lOO-RlOOO-P 



500 



2000 



1000 1500 
Radius ikpcl 
Episodic Jets, Inject. 1000 Myr 




500 



2000 



1000 1500 
Radius [kpc] 
Episodic Jets, Inject. 3000 Myr 




500 1000 1500 
Radius [kpc] 



Figure 4. Time evolution of the logarithm of the mean relative metallicities showing the differences between Runs with a flat (black) 
and a Plumer (red) metal distributions. The dynamical evolution of the ICM metallicity is independent of the metal distribution details 
at the central dominant galaxy. 



runs with the same parameters except for the galaxy metal 
distribution. It is important to note that the total AGN en- 
ergy input to the ICM is the same in all cases. 

Figure [5] shows that although powerful AGN jets are 
transitory phenomena, they produce important dynamical 
and chemical effects on the ICM on megaparsec scales for 
a few gigayears. The large-scale dynamical evolution of the 
ICM and its metals is similar for all our simulations (see 
Figures (2] [3] and |4ll . The initial metal distribution of the 
central galaxy is rapidly lost and the very steep and negative 
gradient of the metallicity rise slowly. The central region, 
r < 250 kpc, remains quite steep specially if metal injection 
is present. From 2 Gyr onwards regions with flat gradients 
develop for r > 500 kpc and for some central regions as well 
(see Figures |3] and [5} . 

The number of jet episodes has an effect from about 
2 Gyr onwards, mainly within the regions r G (200,750) kpc 
and r £ (1500,2150) kpc. The metallicity profiles in Figure [2] 
show that after evolving for 3 Gyr, the ICM metals advected 
by single outburst jets (single jets) reach radii ~ 1800 kpc 
from the cluster's core. However, if intermittent jets are 
present, metals reach ~ 300 kpc further. On the other hand, 
metallicity plots, more signiflcantly from systems with in- 
termittent jets, show reductions in metallicity in the re- 
gion r G (200,750) kpc of an order of magnitude (see Fig- 



ures [2] and O. Similar depressions were found bvl Heath et al.l 
( 2007) in their simulations but it is difficult to find any traces 
of such metal distribution gradients with current X-ray spec- 
troscopic observations of clusters. In our simulations, gas 
and metals are continuously injected into the central domi- 
nant galaxy (see equation ^) and thus the AGN relic bub- 
bles inflated by episodic jets carry with them more metals 
and gas than bubbles formed by single jets do. The ICM gas 
in the surroundings will occupy the space left by the bub- 
bles when they rise, producing central inflows in the cluster. 
These inflows are larger for episodic jets. The combined ef- 
fects of these inflows (replenish central regions with metal- 
poor gas) and the rising bubbles (moving hot gas and metals 
away form the core) result in an intermediate intermedi- 
ate region in the cluster, with low metallicity. Furthermore, 
during the intermittent jets' active phases every outburst 
drives an expanding shock wave leaving turbulent flows in 
its wake. Shocks of secondary outbursts sweepout turbulent 
gas formed by previous outbursts. Thus, the hydrodynami- 
cal structures produced are more complex when intermittent 
jets are present, in agreement with the radi o observations of 
intermittent jets of iBrocksopp et al] (|2007l ). These physical 
conditions facilitate the mixing of gases with different metal 
concentrations, resulting in smoother metallicity profiles as 
shown in the last row of Figures [2] |3] and [5] 
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No Injection 2000 Myr 
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R/R180 



Single Jets 1 000 Myr 




0.2 
R/R180 
Single Jets 2000 Myr 




0.2 0.3 
R/R180 
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0.2 0.3 
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Episodic Jets 1 000 Myr 




0.2 
R/R180 

Episodic Jets 2000 Myr 




0.2 0.3 
R/R180 

Episodic Jets 3000 Myr 




0.2 0.3 
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Figure 5. General evolution of the mean relative metallici ties of our simulations overlaid on the ICM radial metallicity profiles from the 
data composition (Figure 5) bv lLeccardi fc Molendil |200^. 



The metal injection has efltects on the general metal 
abundances present in the ICM, specially at the core 
(see 2M, 2R, 3M, 3R 5M and 5R). The faster the metal 
replenishment time, tc, the higher the metal abundance. Fig- 
ureOshows that the metallicity profiles at radii <75kpc are 
dominated by the central galaxy. On the other hand, an im- 
portant metallicity drop is evident at the core of the cluster 
between 600 Myr and 1 Gyr, for runs with metal injection 
(see 2M, 2R, 3M and 3R). This loss and redistribution of 
metals takes place the moment the equatorial part of the 
jets' cocoon collapses, and the formation of the buoyant 
bubbles of hot gas. In this process some small proportion 
of metals are lost in our simulations through the inner com- 
putational boundary - while others, are lifted to larger radii 
to form the bubbles. However, by comparing Runs with and 
without injection in Figures [5] and [31 from 1 Gyr onwards, it 
seems that the rapid central metallicity changes mentioned 
above (between 600 Myr and 1 Gyr) have no important ef- 
fects on the general metallicity evolution of the ICM. 



The metal distribution details of the central galaxy (ei- 
ther flat or Plumer, see Section 12. 2p have no effect on the 
dynamical evolution of the ICM metallicity on large scales, 
as shown in Figure ID The discrepancies in 4M at 3 Gyr (and 
2M as well) are only due to numerical advection. 



3.2 Comparison with Radio and X-Ray 
observations 

As discussed above, intermittent jets drive turbulence and 
mixing in the ICM. Direct evidence of jets with three out- 
burst episodes was first found i n the FR II radio galaxy 
B0925+420 by iBrocksopp et al. I l|2007l ). The ages of the 
lobes of this radio source (~a few Myr) are comparable to 
the intermittancy timescales we use in our simulations. Fur- 
thermore, there is observatio nal evidence of turbu lence and 
diffusion in the ICM, e.g., see lRebusco et al. I I2OO6I and refer- 
ences there in. However, a quantitative study of the "nested 
cocoons" and the ICM turbulence and diffusion is beyond 
the scope of this letter. 

In Figure [S] we compare the m ean metallicity pro- 
files fr om the X-Ray data presented in lLeccardi fc Molendil 
(|2008h and references there in, to our simulations at differ- 
ent epochs. It is evident that for radii < 0.12 Rigo our sim- 
ulations present pronounced metallicity gradients that do 
not follow the general trends present at larger radii. In our 
simulations this inner region is dominated by the metals of 
the central galaxy (see Figures [2] and [3]) . We do not expect 
our simulations to reproduce the metallicity gradients on 
small scales, where the profiles are likely to be shaped by 
less-energetic phenomena than AGN activity such as galac- 
tic winds. Away from the core however some of our simu- 
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Table 2. Ranges of agreement between our simulations and the 
observations from lLeccardi fc MolendU l|2008l ) . at 3Gyr. 



Simulation 


Approximate radial 


name 


range [Riso] 


lO-RInf 


(0.18,0.45) 


lO-RlOOO 


(0.26,0.40) 


I0-R500 


(0.32,0.45) 


IlO-RInf 


(0.27,0.42) 


1100- RInf 


(0.31 ,0.45) 


IIOO-RIOOO 


(0.29,0.45) 


1100- R500 


(0.25,0.45) 



lations develop metallicity gradients comparable, over radii 
from 0.1 to 0.4R/Ri8o, to those in the data (see Figure [5}. 
Table 2 summarises the radial ranges over which the mean 
metallicity of our simulations is in agreement with experi- 
mental data. The best fit is obtained with Run lO-RInf, with 
a single jet outburst and no metal injection. 



4 CONCLUSIONS 

We present 2D axisymmetric hydrodynamical simulations of 
the metal advection in the ICM, driven by the jets of a radio- 
loud AGN. We implement the ICM with a cooling flow, and 
the central cD galaxy with a metal distribution and contin- 
uous metal formation, and powerful intermittent jets. The 
jets produce the early and late features typical of canonical 
models, simulations and observations of F R II objects (e.g. , 
see[C hurazov ct al. 20 01; Alexander 2006: . Heath et al.ll2007l : 
iBrocksopp et al.il2007l ). We model the metals using tracer 
fields and explore the consequences of varying the jet inter- 
mittency, the galaxy metal distribution profile and the metal 
formation rates. We have found the following. 

Although powerful AGN jets are transitory phenomena 
they produce important dynamical and chemical effects on 
the ICM on megaparsec scales for a few gigayears. 

The convective flows formed long after the AGN active 
phases drag ~95% of the initial metals in the central galaxy 
to distances > 1.5 Mpc. 

We implement intermittent jets with the same total 
energy input to the ICM. The general metallicity evolu- 
tion for all our simulations is very similar on large-scales: 
after 2 Gyr, approximately flat metallicity gradients are 
formed from radii ~ 700 to ^ 1300 kpc. Metallicity drops are 
formed (found by iHeath et al.l l|2007l ) as well) of factors 
of ~ 15 at r £ (200,750) kpc. These features might be corre- 
lated to the mass of the gas and metal advected by the relic 
bubbles, but high resolution X-ray observations are needed 
to see traces of such metallicity gradients. Furthermore, in- 
termittent jets distribute metals ~ 300 kpc further than sin- 
gle outburst jets do, and also facilitate the mixing of gases 
with different metal concentrations and induce turbulence 
in the ICM. 

The metal injection has effects on the general metal 
abundances present in the ICM: the faster the metal re- 
plenishment time the higher the metal abundance. For 
radii < 0.12 Rigo our simulations are dominated by the met- 
als of the central galaxy and do not follow the general trends 
present at larger radii on either our simulations or the X- 



Ray observations. Since metal production and jets are real- 
istic in active galaxies, metalliciy gradients in the very cen- 
tral regions of cool core clusters are likely to be shaped by 
less-energetic phenomena. On the other hand, away from the 
core, some of our simulations (specially the one with a single 
jet outburst and no metal injection) develop metallicity gra- 
dients in agree ment with, over particular ra dial ranges, the 
observations of iLeccardi fc Molendil l|2008l ) and references 
there in. 

We model massive central galaxies with the same total 
metal mass in all cases. The metal distribution details of the 
central galaxy have no effect on the dynamical evolution of 
the ICM metals. 
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